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I present a review of selected topics in the computation of heavy flavour cross sections in photon-hadron 
collisions. 



1. Computation of cross sections at fixed 
order in QCD 

A heavy quark Q is by definition a quark whose 
mass mq is much larger than Aq^u. Since iuq 
is different from zero, it cuts off the singulari- 
ties due to the collinear emission of gluons from 
the heavy quark lines, thus allowing the defi- 
nition of open-heavy-quark cross sections, that 
are infrared-finite order by order in perturbative 
QCD. The hard scale of the production process 
is of the order of tuq. Therefore, one expects 
the perturbative predictions to be very reliable in 
the case of top, under control in the case of bot- 
tom, and affected by large uncertainties in the 
case of charm. This is more so, since the smaller 
the mass, the bigger the impact that effects of 
non-perturbative origin will have on the predic- 
tions of physical observables. The cross section 
for the production of QQ pairs in on-shell photon- 
hadron reactions is usually written as the sum of 
two terms: 

denoted as pointlike and hadronic components re- 
spectively. Neither of the two is well defined be- 
yond LO in perturbation theory. This fact is for- 
mally taken into account in eq. by the depen- 
dence upon the scale fi^, which is related to the 
subtraction of singularities due to the collinear 
splitting of the incoming photon. The pointlike 
and hadronic components are not constant with 
respect to a variation of fi^ , while their sum is (up 
to higher order terms in perturbation theory) . On 
the other hand, the pointlikc and hadronic com- 



ponents are separately constant with respect to a 
variation of all the other mass scales which enter 
eq. (Q). The contribution of the hadronic com- 
ponent will be larger (relative to the pointlike 
component) the smaller ttiq. However, it has to 
be remarked that kinematical configurations exist 
where the hadronic component is basically negli- 
gible: typical examples are the high-p,, region of 
the emitted heavy quarks, and small center-of- 
mass energy collisions (like those at fixed-target 
experiments). This is due to the fact that photon 
densities f^^^ (x) are vanishing faster than hadron 
densities for x —> I. 

The task for theorists is to calculate the par- 
tonic cross sections a^i and aij at the highest 
possible order in perturbation theory. In ref. [ [^] 
(for the pointlike component) and in ref. [ ||] (for 
the hadronic component) results accurate to NLO 



(ae 



and ttg respectively) have been obtained. 



At that time, predictions were only available for 
total cross sections and single-inclusive spectra. 
Later, computations were performed which 
allowed the definition of cross sections fully exclu- 
sive in terms of the variables of Q, Q, and of the 
accompanying jet. The matrix elements relevant 
for QQ production are mostly convenient calcu- 
lated in a modified MS scheme [ ^ , in which the 
renormalized as is defined through the following 
equation (up to terms of order a^): 
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^Talk given at DIS99, 19-23 April 1999, Zeuthen, DE. 



where Nif is the number of flavours lighter than 
Q (i.e., 3 for charm, 4 for bottom and 5 for top). 
The first term in square brackets in eq. (|) is the 
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same that appears in the standard MS scheme; 
thus, hght flavours are dealt with like in stan- 
dard MS. The second term in square brackets is 
on the other hand peculiar of heavy flavour pro- 
duction, and is designed in such a way that, for 
momenta much smaller than TOq, it cancels ex- 
actly the contributions of heavy quark divergent 
loops. Since this second term is only expressed 
in terms of a^'"''^ , it follows that the renormaliza- 
tion group equation one derives is equal to that of 
the standard MS with Ni f flavours. This is con- 
sistent with the fact that no contribution to the 
cross section comes from diagrams where a heavy 
quark is present in the initial state. 

The cross sections computed in this way are ex- 
pected to give sensible predictions for momenta of 
the order of ttiq. At the borders of the phase 
space, or for high energies, some logarithmic 
terms may grow large and spoil the convergence 
of the perturbative series. Example of large loga- 
rithms are: log(5'/mQ) (small-x effects: they are 
important when the center-of-mass energy of the 
colliding particles is large); loglpT^^ /itiq) (due to 
the emission of soft gluons, they are important 
at the edges of the phase space or at the thresh- 
old); log{pT/mQ) (due to the emission of coUinear 
gluons from the heavy quark line, they are impor- 
tant when Pj, ^ mo). It has been shown that the 
impact of small-x effects for charm and bottom 
production at present colliders is moderate. Also, 
soft gluon resummation gives small enhancements 
with respect to NLO results for total cross sec- 
tions in the case of top production at the Teva- 
tron and of bottom production at HERA-B. In 
the following, I will concentrate on coUinear log- 
arithms, which are phenomenologically relevant 
for charm photoproduction at HERA. 

2. The large transverse momentum region 

If Pt ^ , mass terms in the matrix elements 
are suppressed by powers of wlqIpt. Therefore, 
when computing matrix elements one is entitled 
to treat the heavy quark as massless [ H, ^ and 
use the standard MS scheme (this is the reason 
why this kind of calculations are often - incor- 
rectly - referred to as "massless"). This implies 
that it is mandatory to consider diagrams where 



the heavy quark is also present in the initial state. 
Furthermore, the flnal state Qg coUinear singu- 
larities are no longer cut off by the quark mass. 
The subtraction of these additional singularities 
is performed within the framework of the factor- 
ization theorem; one writes the cross section for 
the production of a heavy quark (in any kind of 
hard collisions) as follows: 

daQ^d<j,{ti)®D'^\ii), (3) 

where i can be a gluon, a light quark or the heavy 
quark Q, and dai is the cross section for the pro- 
duction of the flavour i. Dq are the so-called per- 
turbative fragmentation functions (PFF), whose 
bare form contain a divergent term which cancels 
the final-state coUinear singularities mentioned 
above. As far as the ^ dependence is concerned, 
the PFFs obey the usual Altarelli-Parisi equa- 
tions. However, at variance with the case of par- 
ton densities, the PFFs are fully calculable in per- 
turbation theory, for scales ^ = /ig of the order 
of the heavy quark mass (this prevents Dq\ij,q) 
from containing large logarithms). Such a calcu- 
lation has been performed in ref. [ |^ , up to terms 
of order a|. At this order, d'^q'' and Dq^ are 
non- vanishing, while d\IP is zero. The main fea- 
ture of eq. (^) is the following: the scale is of 
the order of Pt, which is the hard scale of the 
process. Since in the computation of matrix ele- 
ments the mass of the heavy quark has been set to 
zero, d<Ti does not contain any large logarithm of 
PtI^q- All the mass effects are included in ZJq''; 
however, here the large logarithms are resummed 
by the Altarelli-Parisi equations, in the evolution 
from the scale /xq to the scale ^. 

If the data are presented in terms of heavy- 
flavoured hadrons Hq, theoretical predictions 
must be given for such quantities; these are 
calculated starting from the heavy-quark cross 
section given in eq. (^), and convoluting it 
with a non-perturbative fragmentation function 
(NPFF), that accounts for the hadronization pro- 
cess Q — > Hq-. 

da„^ = dcjQ ® Di?^ (4) 

(of course, the convolution with a NPFF can 
be also performed starting from a fixed-order 



3 



QCD cross section, like the one in eq. (|i|)). The 
NPFFs can not be calculated in perturbation the- 
ory; however, they are supposed to be universal. 
Therefore, they can be determined by a fit to ex- 
perimental data - typically, from e^e~ collisions 
- and eventually used to obtain theoretical pre- 
dictions for various production processes. Recent 
fits are presented in refs. [ ^, ^; it turns out that 
the NPFF is much harder than in previous deter- 
minations based on LO QCD analysis. A different 
procedure has also been proposed [ |^. Instead 
of considering PFFs and NPFFs, one can intro- 
duce a fragmentation function which describes at 

(i) 

the same time the hard and soft physics: Dh^ — 

ZJq'' ® . The idea is then to parametrize D^g , 
fit it to the data, and then evolve it to the appro- 
priate scale with Altarelli-Parisi equations. How- 
ever, since the NPFFs considered in practice do 
not depend upon a mass scale, there is no differ- 
ence in evolving only the PFF or the convolution 
of the PFF and the NPFF (this can be trivially 
shown by writing the Altarelli-Parisi equations in 
Mellin space). Therefore, ref. [|| and on/?/ dif- 
fer in the choice of the initial conditions for the 
evolution of the fragmentation function. While in 
ref. [ D the as terms of Dq\iiq) are taken into 
account, they are not included in ref. [ ^. Thus, 
ref. [ P only includes mass effects through a fitting 
to data, and does not exploit all the perturbative 
results available. 

3. Matching large and small momentum 
results 

As discussed in the previous sections, if is of 
the order of mg, then fixed-order QCD compu- 
tations have to be used; on the other hand, if 
is much larger than mg, resummed calculations 
are mandatory. What is missing, unfortunately, 
is a quantitative definition of what "much larger" 
means. HERA data are compared to resummed 
calculations for Pt{D*) down to 2-^3 GeV, which 
means Pt (c) on average smaller than 4 GeV (here 
and in what follows, I take Pt{D*) — {z)pt{c), 
where {z) is the average z obtained using a Peter- 
son fragmentation function with e — 0.03). The 
ratio 4 GeV/mc does not seem to be large enough. 



Any agreement between data and resummed cross 
sections in this Pt region should therefore be re- 
garded as incidental. However, it is clear that 
there is an intermediate region in p^ where both 
fixed-order and resummed calculations are almost 
equally reliable (or, if one adopts a pessimistic 
point of view, neither of the two is meaningful). 
Therefore, it is conceivable to have a matched cal- 
culation, which returns a fixed-order result at low 
Pt , a resummed result at large Pt , and smoothly 
connects the two approaches in the intermediate 
region. Unfortunately, such a calculation does 
not exist for photoproduction. However, results 
are available for hadroproduction [ [l^, [ill . It is 
instructive to see how this matter is dealt with 
in ref. [ since there is no conceptual differ- 
ence with respect to the case of photoproduction 
(to recover the latter, simply substitute a" with 
Q!emQ!"~^ in the following equations). The fixed 
(NL)-order and (NLL) resummed cross sections 
are written as (/i is a scale of the order oi p^): 




+ o{al[as\og^jy (6) 

The coefficients A and B contain logarithms of 
Pt/^q- Eq- (^ is defined up to power-suppressed 
terms, not explicitly indicated. We can now de- 
fine the matched cross section: 




where the function G is arbitrary to a large ex- 
tent, but must approach 1 when mq/pT 0. 
If Pj. ~ ttiq , the quantity in curly brackets van- 
ishes, and the matched cross section approaches 
eq. (^. If Pt ^ mq, the matched cross section 
approaches eq. (|^), since logarithmic terms in A 
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and B dominate (by construction, these terms are 
exactly those missing in the sums of eq. (|^) with 
respect to the sums of eq. (^)). Notice that this 
is true regardless of the form of G, provided that 
it satisfies the condition given above. One could 
therefore choose G = 1. However, as shown in 
ref. [ 0, this choice can lead to instabilities in 
the numerical computations. In ref. the fol- 
lowing form has been adopted: 



G{mQ,pT) 



(8) 



where /3 is a suitable number. From numerical 
studies performed for bottom production at the 
Tevatron, it turns out that a sensible choice is 
/? = 5. This means that, in the matched cross sec- 
tion, the whole tower of logarithms (except those 
with coefficients ag, ai and Bq) is suppressed by 
a factor 0.5 for — ^ttiq, and is almost fully 
included (G = 0.9) for = ISmg. If we as- 
sume that these conclusions can be also applied 
to charm production at HERA, we conservatively 
conclude that a resummed cross section can be 
safely compared to the data starting from a Pt (c) 
of the order of 22 GeV (here, rric — 1.5 GeV), 
that is for PriD*) of the order of 16 GeV. Be- 
ing less conservative, we can say that the lowest 
Pt{c) value for which a comparison between the 
data and a resummed cross section is meaningful 
is of the order of 8 GeV {pt{D*) of the order of 
5^6 GeV). 

4. Conclusions 

I presented a review of results relevant for 
the production of heavy flavours in high-energy 
photon-hadron collisions. At low and moderate 
transverse momenta, where the bulk of the total 
cross section comes from, fixed-order QCD com- 
putations (available to NLO accuracy) are appro- 
priate. HERA data on charm production have 
renewed interest in the resummation of large log- 
arithms of the transverse momentum. A brief 
discussion on the techniques designed to perform 
this resummation has been given. Finally, it has 
been argued that most of the currently available 
data for D* production at HERA have too small 
a Pt to be safely compared with a resummed cal- 



culation. A definite statement on this issue will 
however only come when a matched calculation 
will be available, which combines the results of 
fixed-order and resummed computations. 
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